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Abstract 



The time spent by an interacting Brownian molecule inside a bounded 
microdomain has many applications in cellular biology, because the 
number of bounds is a quantitative signal, which can initiate a cas- 
(— I , cade of chemical reactions and thus has physiological consequences. In 
Qh| the present article, we propose to estimate the mean time spent by a 
Brownian molecule inside a microdomain which contains small holes 
on the boundary and agonist molecules located inside. We found that 
the mean time depends on several parameters such as the backward 
binding rate (with the agonist molecules), the mean escape time from 

■ the microdomain and the mean time a molecule reaches the binding 
, sites (forward binding rate). In addition, we estimate the mean and the 

\0 variance of the number of bounds made by a molecule before it exits 

These estimates rely on a boundary layer analysis of a conditional mean 
^ first passage time, solution of a singular partial differential equation. In 

■ particular, we apply the present results to obtain an estimate of the 
mean time spent (Dwell time) by a Brownian receptor inside a synaptic 
domain, when it moves freely by lateral diffusion on the surface of a 



neuron and interacts locally with scaffolding molecules. 

Introduction 



Biochemical reactions in cellular micro domains involve, in general, a small 
number of molecules that can bind to agonist molecules, confined in some 
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subregions. A microdomain Q is defined as a bounded domain where a large 
fraction of the boundary dQr is reflective and a small part dQa of it is absorp- 
tive, which allows molecules to enter or/and exit. Under normal physiological 
conditions, a molecule can be trapped inside a microdomain for a period long 
enough, compared to other time scales such as free diffusion or binding time. 
Because a molecule can be trapped for a long time, many chemical bounds 
can be made, before it exits. From a physiological stand point, an interesting 
property of such microdomain is that depending on the number of chemical 
bounds made by a molecule, a cascade of chemical reactions can be initiated, 
which ultimately affects some physiological properties. Our interest here is 
to estimate the number of chemical bounds made by a Brownian interacting 
molecule inside fl. This number depends on the geometry of the microdomain 
and the distribution of agonist molecules. We also provide an estimation of the 
mean time spent (Dwell time) by a molecule inside the microdomain, denoted 
by E{td)- For that purpose, we derive an asymptotic formula for the Dwell 
time, the mean and the variance of the number of bounds made before the 
molecule exits, as the ratio tends to zero. Finally, we derive a formula for 
the Dwell time of a molecule inside fi, when a steady state flux of molecules 
is maintained fixed at the absorbing boundary. The role of the flux is to fix 
the number of free binding sites. Using the present method, we also obtain 
an estimate of the forward binding rate constant. Historically, the theory of 
chemical reactions at a molecular level limited by diffusion has been developed 
by many authors, to quote but a few [1, 2, 3, 4, 5]: using the classical theory of 
diffusion and interactions with binding sites, various rate constants were com- 
puted. Recently, using averaged equations, chemical reactions in microdomains 
have been described in [6, 7]. Chemical reactions in closed microdomains were 
studied in [8], where we obtain, in particular, some estimates on the mean of 
the variance of the number of bound molecules in a steady state regime. In 
[9] , an asymptotic estimate of the mean time it takes for a Brownian molecule 
to escape an empty domain through small openings located on the boundary, 
was obtained using a new type of singular perturbation problem. More specif- 
ically, if D denotes the diffusion constant, \Vt\ the volume of the domain Vl and 
£ — -i^^ << 1 is the ratio of the absorbing to the total boundary, then for 
e small, the leading order term of the mean escape time t{x) (for a molecule 
starting at position cc, far from the entrance) is given by 



In the first approximation, the mean time t{x) does not depend on the initial 
position X and will be denoted by r. In this article, we obtain an explicit 
asymptotic estimation of the Dwell time Eiru) as a function of the character- 
istic sizes of the domain Q, the size of the small openings dO^ai the number of 
the binding molecules. More specifically, E{td) is given by expression (2.12), 
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which depends on the mean time (r) to exit when no binding occurs, the mean 
time (T) to enter into the binding site area, the mean probabihty to bind 
before exit and the backward binding rate k_i. We get 



It is well known from the theory of chemical reactions that the backward 
binding rate A;_i depends only on the local interactions between two interacting 
molecules. If denotes the activation barrier, kT^. is the energy due to the 
temperature, the Arrhenius law states that: 



where C is a constant that depends on the temperature Te, the electrostatic 
potential barrier A£' generated by the binding molecule and the friction coef- 
ficient [10]. In the first part of the paper, we derive equation (1.2) by count- 
ing the number of bounds between the Brownian molecule and the agonist 
molecules, before the Brownian molecule exits the domain. In the second 
part, we derive some asymptotic estimates of the quantities (T), (r) and rris 
as a function of the geometry, when the radius 5 of the binding site tends to 
zero. Although the present computations are carried out in two dimensions, 
they can be extended to dimension 3 by using the techniques developed in [11]. 
Finally in the last part, we apply the present computations to study chemical 
reactions occurring in synaptic microdomains: the chemical reactions arc the 
binding of receptors with the scaffolding molecules. It is indeed of great interest 
to analyze the mechanism that regulates the number and the type of receptors 
at a synapse, because receptors control the synaptic weight. Any fluctuations 
of the number results in a variation of the synaptic weight and affects the 
reliability of the synaptic transmission. Moreover, certain experimental proto- 
cols have lead to a Long Term Potentiation of a synapse, a mechanism which 
is associated with a change of the number and the type of certain receptors 
[12, 13]. The regulation of synaptic plasticity is a fundamental process under- 
lying learning and memory [12, 13] and recently, single molecule tracking has 
revealed that the number of postsynaptic receptors, which participate to the 
synaptic transmission, is not fixed but it changes due to constant trafiick of 
receptors on the surface of neurons. Receptors move in and out from synaptic 
regions [14] [15] and following these observations, many questions have been 
raised: in particular, what determines the time spent by a receptor inside a 
synapse? How receptors can be stabilized inside a synapse? How long they 
stay inside synaptic microdomains? Such questions are partially answered in 
the present paper. In particular, our computation of the Dwell time of a re- 
ceptor inside a specific microdomain, called the Postsynaptic density (PSD) 
takes into account the interaction with the scaffolding molecules, which was 
ignored in [9]. 




(1.2) 



/c_i = Ce , 



(1.3) 
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1.1 Molecular dynamics in a microdomain 



The dynamics of a molecule or a protein moving on the surface of a cell is 
usually described in the Smoluchowsky limit (large friction) of the Langevin 
equation [10]: for a molecule of mass m, described by its position X at time 
t. with a friction coefficient 7, moving inside a potential well V, the Smolu- 
chowsky limit of the Langevin equation is 

7X + Vy(X) = ^2^ti;, (1.4) 

where = — and w is a Gaussian random variable with variance 1 and mean 
0. In a microdomain Q, where a large fraction of the boundary is reflective dflr 
and a small part of it is absorptive dfla, the probability density function (pdf) 
p to find X at time t in the surface element x + dx satisfies the Fokker- Planck 
Equation (FPE) 



dp{x, t) 



= DAp{x,t) -W ■ {WV{x)p{x,t)) ioi X en (1.5) 



dt 

J{x, t)-n ^ for a; e dflr (1-6) 

p{x,t) = Qioixedna (1-7) 

where D — ^e^ is the diffusion constant, n is the external normal at the 
boundary, the flux J is given by 

J{x, t) = -DVp{x, t) + VV{x)p{x, t). (1.8) 

We denote by the first time the molecule arrives at the absorbing boundary 
dQa, when it started at position x. It is well known [10] that the mean first 
passage is the expectation of the time and is given by 

/"OO J /-oo 

E^if") = J t—Pr{t''<t}dt^j Pr{t''>t}dt 
p{y,t\x)dydt, 



where p{y, t\x) is the pdf of the process X, conditioned on the initial position 
X, that is, as t goes to zero, 

p{y,t\x)^5{x-y), (1.9) 

where 6 is the Delta-Dirac function. In equation (1.5), V represents the po- 
tential wells generated by the binding molecules inside the domain Q. It is, in 
fact, the sum of the local potential wells, supported inside a ball of finite radius 
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generated by the binding molecules. Usually the binding molecules are scat- 
tered inside the domain Q, but in the present model we replace the scattered 
distribution of binding molecules by a simplified distribution, where we imag- 
ine that all the binding molecules are located inside a compartment D{6) in Q. 
In the present description, the potential V becomes an effective potential de- 
fined in D{S), whose characteristics should be such that the classical chemical 
reaction theory is recovered. More specifically, we can define the microdomain 
Q containing the binding domain D(d), which replaces the domain with many 
scattered binding sites: this simplified domain made of two compartments is 
called the homogenized microdomain and is described as (see figure 1) 

1. An internal compartment, which is described as a disk D[S) of radius S. 
This disk represents the region where the binding sites are located. In- 
stead of using the dynamics associated with equation (1.4), we describe 
the entrance and the exit of a molecule inside D{d) using a Poissonnian 
description, where the mean can be related to the properties of the po- 
tential well. For that purpose, we recall that a chemical reaction with 
a binding molecule is described as the arrival of a Brownian molecule 
inside the disk D{6). The release process is modeled as the escape of 
the molecule from the potential well V and is described by the chemical 
reaction 

h 

R + S # RS (1.10) 

where y~[ is the mean binding time and depends only on the local po- 
tential well, generated by the binding molecules [16] [10]. 

2. An external compartment separated from the rest of the biological envi- 
ronment by the boundary dQr, containing small openings dQa- Molecules 
can enter or exit through the openings and thus can be exchanged with 
the rest of the cellular medium (see figure 1). The dynamics of a molecule 
in that compartment is described as pure Brownian until it escapes. 



1.2 Time spent by a molecule inside a disk containing 
a small opening 

We estimate asymptotically the Dwell time E{td) of a molecule inside a domain 
Q, when the domain Q is a disk D{R) of radius R. The case of a general domain 
is left open. By definition, E(t£)) is the mean time to exit, averaged over an 
initial uniform distribution. We make the following assumptions: the ratio e of 
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Figure 1: Model of a cellular microdomain in two dimensions. The domain is 
a disk D{R) of radius R, made of two compartments: an inner disk D{S) of 
radius 5 and the annulus As = D{R) — D{d). A molecule moves by Brownian 
motion inside As until it hits D{6) or the absorbing boundary dQa- When the 
Brownian molecule enters into D[6), which represents the domain of chemical 
reactions, while it reacts with an effective binding molecule during a mean time 
(the inverse of the backward binding rate) its movement in D(S) is frozen. The 
molecule is then released uniformly inside the annulus As = D{R) — D{5). This 
scenario repeats until the molecule hits the absorbing boundary dVta, where 
the molecule is finally removed. 

the absorbing to the reflective boundary of D[R) is small << 1. The absorbing 
(resp. reflecting) part is denoted dVla (resp. dVLr)- The computation of the 
Dwell time will be made using the homogenization version of the domain, 
described in the previous paragraph: the binding molecules are located in the 
small concentric disk D(5), where S << R. Each time a molecule enters into 
D((5), it can be bound with a binding molecule and will be released in the 
annulus As(R) = D{R) — D{6). Since R is a flxed parameter, we will denote 
the annulus by As. We assume that the boundary dD{5) is always absorbing. 

2 The Dwell Time E{td) 

The goal of this section is to derive a general expression of the Dwell time 
E{tb)- The movement of a molecule is described here as Brownian (see [17, 14] 
for the case of a synapse) with a diffusion constant D and X* {i) is the position 
of a molecule at time t starting at position x. Eijr)) is the mean time a 
molecule remains inside the domain Vt before it exits, averaged over a uniform 
initial distribution inside the annulus As = D{R)—D{6). Inside As, a molecule 
diffuses freely until it hits dD{5). Inside D{5), a bound molecule is kept flxed 
for a mean time before being released inside As at a random position, 
distributed uniformly. To estimate the Dwell time E{td), we flrst consider the 
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Dwell time E(t^), conditioned on the initial position located inside As. We 

do not consider here the entrance of the Brownian receptor inside the domain, 
which would require to use a Langevin formulation of the dynamics [10]. We 
derive a formula for E{t^) by counting the number of times a molecule enters 
into D{S) before it exits: either a molecule exits with no bindings or many 
bindings occur, before the receptor leaves D{R). Let denote the first 
passage time the trajectory X^{t) hits the absorbing boundary dQa- Similarly 
we define by Tj the first passage time of the trajectory X^{t) to the boundary 
of D{6). We define now the probability ps{x) that a trajectory leaves 
the domain Ag before any bounds occur, explicitly 

ps{x) = Pr{r^ < r|}. (2.1) 



2.1 The general formula for the mean Dwell time E(td) 

To estimate E{td), we consider the different random events related to the 
number of times the Brownian molecule enters into the region As, before it 
exits. is the event that a receptor exits the microdomain As when no 
bindings occur, starting at position x. F^ is the event that a receptor starting 
at X, enters the domain D{5) only once and then exits without entering again 
in D{5). Similarly, we define the event F^ that a receptor starting at position 
X exits after exactly n bounds. The probability that the molecule has left the 
domain before time t is Pr^T'"" < i}. The mean time is given by 

E{t^) = / t^Pr{T^<t}dt (2.2) 

Jo dt 

= tf^Pr{r'^<t,F^}dt. (2.3) 

We use Bayes formula to get 

E{t^) = E{r^\F^) Pr(Fo^) + i;(T^|Ff ) Pr(Ff ) + E{r^\F^) Pr(Ff ) + ...(2.4) 
where 

Pr{T^ < t\F^} is the probability that < t conditioned to the event that 
exactly n bounds are made before the molecule exits the domain D{R). When 
a molecule detaches from the domain D{6), in this model, it is released inside 
the domain As according to a density distribution p{x): immediately after the 
mean duration 7^, the molecule is released inside the volume element x + dx 

K—l ' 
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with a probability p{x)dx. To compute the first term in equation (2.4), we 
notice that 

E{t^\F^) Pr(F-«) = E{T^\T^ < Ti)ps{x). (2.6) 

The second term E{t'''^\F^) Pt{F^) is the sum of two terms. First the mean 
time it takes for a molecule to reach the boundary dD[S) and stays there for 
a mean time ^ and second the mean time to go from a point Xi, where 
the molecule starts afresh, to the absorbing boundary. If we note that the 
probability to reach dD{5) before dQa starting from cc is 1 — ps{x), then we 
have: 

E(r^|Ff )Pr(Ff) = {I - ps{x)) x 
(^{T||T| < T^} + -^) / p{x,)ps{x,)dx,+ I E{T^^\T^^<Tijp{x,)ps{xi)dxi 

K-l J As J As 

Similarly, the general term where exactly k bounds are formed is computed, 
with the notation 

t^{x) = E{Ti\T^<Ti} (2.7) 
t'(x) = E{Ti\Ti,<T^}. (2.8) 

We get, 

,A 



J As Jas fc-i 

n-l 

JJ(1 - psiXi))p{Xi){l - ps{x))ps{Xn)p{Xn)dXi . . . dx^- 

i=l 

= {l-ps{x))\^t^{x){l-ms)''-^ms 

+ (n - 1)(1 - msY'-^ms / f{y){p{y) - P5{y)p{y))dy 

Jas 

+ ^(1 - msT-'ms + (1 - msT-' t^{y)ps{y)p{y)dy} 



where 



ms= ps{y)p{y)dy. 

JAs 



'As 

Formula (2.4) is in fact a geometric sum and using the generic expression of 
the series, we get 

E{t^) = t\x)ps{x)-r^-^^^{t'{x)-r I t^{y)ps{y)p{y)dy (2.9) 

^i<5 JAs 

+ —if t'{y){p{y)-ps{y)p{y))dy + ^)\. 
ms Jas k-i ) 



When p is uniformly distributed, formula (2.9) can be simplified. Using that 



ms = 



As 



1 




ps{y)dy, 



where 1^451 is the volume of Ag and if we define the two mean times 



(r) 




t^{x)ps{x)dx 



(2.10) 



lAsPsiy)dy 



and 



(T) 



(2.11) 



then the Dwell time E{td) is given by the average of E^t""") with respect to 
the uniform distribution and we get 



In the next paragraphs we obtain an explicit asymptotic estimate of for com- 
ponent of formula (2.12). 

Remark. It seems artificial to free a bounded molecule not immediately 
at the boundary of the domain Dg. But in fact, the release process can be 
modeled by re- initiating the molecule at the boundary of Dg, otherwise it would 
immediately return to a bounded state. To avoid this unrealistic behavior, we 
have used the distribution function p{y) to model the releasing process. A 
more accurate scenario would require to use a Langcvin description, which 
accounts for the acceleration of the molecule: when a molecule is released with 
a random Gaussian initial velocity, it can travel up to a certain distance, before 
its velocity reaches the value of mean thermal velocity of free molecules. 

2.2 The mean number of bounds before exit 

To estimate the mean number of bounds made by a molecule before it exits, 
we use the analysis of the previous paragraph. The probability that a molecule 
starting at x does not bind before exit is exactly Pr(FQ^), while the probability 
that exactly k bounds are made is Pr(F^). We define the average probabifity 




(2.12) 



by 
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The probability Pt{F^) can be expressed in terms of the conditional proba- 
bility p as follows 

Pr(Ff)= / ps{y)p{y)dy{l-ps{x)) 

J As 

and the average probability Pr(Fi) is given by 

Pr(Fi) = / p{x)Fr{F^)dx = / ps{y)p{y)dyp{x){l - ps{x))dx ^ ms{l - ms). 

J As J As J As 

More generally, 

Pr(Ffe)= / ... / psiyk)p{yk)dyk{ps{yk-i)piyk-i)dyk-i)-{'i--psiyi))p{yi)dyi{i-ps{x)) 

JAs JAs 

and after some integrations we get 

Pr(Ffc) = ms{l - ms)''. 

The mean Mb and the variance Vb of the number of bounds before exit are 
given by the well known formula of geometric probability, 

Mb = J]A;Pr(Ffe) = ^A;m5(l-m5)'= 

1-ms ^ lAsi^-Ps{y))dy 
ms J^^ps{y)dy 

Vb = J2^'Fr{F,) - {J2kPr{F,)r = ■ (2.14) 

Later on we will give an asymptotic expansion of nis, Vb and Mb as a function 
of the parameter S. 



3 Estimation of the probability 

ps{x) = Pr{T^ < Tg}. 

To obtain an asymptotic expansion of the probability ps{x) that a molecule 
exits before it enters into the domain D{6), that is, before it binds to any 
binding sites, we use the notations of the previous section. denotes the 
first time a molecule starting from x hits the absorbing boundary. We assume 
that when a molecule exits the domain Q, it does not come back. is the 
first time a molecule hits the inner circle dD{5). The probabihty qs{x) — 
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1 - psix) = Pr{r| < T^} of the event > T§. satisfies an elliptic partial 
differential equation [18], with mixed boundary conditions given by 



Ag5 



■ on As, 
on dVtr, 

qs{x) = on dila, 
qs{x) — 1 on dD{5), 



(3.i; 



where dfta is the small opening located on the external boundary of Q, 
is the remaining part of the external boundary, which is reflecting. In po- 
lar coordinates {r,6), the portion of the boundary dQa is parameterized by 
1^^ — 7r| < e. By using the particular geometry of the annulus As, an explicit 
estimation of the solution can be obtained by using spectral methods, devel- 
oped in the context of mixed boundary value problems [19, 20, 21]. Using 
the method of separation of variables, the solution qg of problem (3.1) has the 
general form 



Qs 



oo (- 

^) = y + E K(^)" + ^°«(^^) + "log ©■ (3-2) 



n=l 



We wish now to estimate the coefficients a„ and 6„. We denote P — Prom 
the boundary conditions on dD{6) and on r = i?, we get 

[any 

n=l 



cos{n9) = 1, 



f.E[-(^)"-^ 

n=l 

^ VR R 

n=l 

oo 

1 + E [«n + cos(ne) - Q;log(/3) 



cos(n^) -I- — 
R 



for 16* - 7r| > £ 
for|e-7r|<£. 



(3.3) 
(3.4) 
(3.5) 



n=l 



Using equation (3.3), we get the following identities: 

ao = 2 

bn = —an/3" for n > 1. 

Using the identities above and (3.4) and (3.5), we obtain the double series 
equations 



n=l 



1 + e: 



n=l 



n % + cos(n^) + I = 0, for |^ - 7r| > £ (3.6) 

an - a„/?2"j cos(n^) - alog{p) = 0, for |^ - 7r| < e. (3.7) 
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Substituting c„ = a„(l + Z?^") and — j^^^ equations (3. 6), (3. 7) have the 
following form 



where 



f + E~ 1 imr cos(n^) = 0, ee[7:,7:-e] 



a + Y,'^=i cos(n6') =0, ^ e [0, tt - , 



| = l-alog(/3). 



(3.8) 



(3.9) 



(3.10) 



The asymptotic solution of equations (3.8)-(3.9) uses the double series expan- 
sion, developed in [19, 20] and used in [9] in the context of a small opening 
asymptotic. In Appendix A, the general solution is given. By using these 
results, we have the following expression for the coefficient Cq 



Co = —2a 



21ogJ + 21og2 + 4/32 + 0(/32,£)]. (3.11) 



Using equation (3.10) and (3.11), we get the expression 

1 



a 



( log(i) + [2 log I + 2 log 2 + 4/32 + 0{P^ £)] ) 



(3.12) 



To remember that a depends on (3 and e, we denote it a{f3,e). For e fixed 
and S small, the other coefficients are estimated by using the expression of c„ 
(given in the appendix) by 



where c„ depends only on n and 



6„ = -c„/3"~0(a(/3, £)/?-). 



(3.13) 

(3.14) 
(3.15) 



These estimates show for S small that the leading term of q is given by : 

l + a((3,e)\og(r/5) + 0((3) for r ~ 5 

1 + a{(3, e) log{r/S) + 0{a) for r ~ i?. 



qs{r,9) 



(3.16) 
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The probability ps{x) that < is given by ^^(cc) = 1 — qs{x) and the 
average over a uniform distribution in formula (3.2) using (3.16) gives 



m = {ps) = 



Ia{S) P5{r,0)rdrde 



= -2 



a{(3, e) \og{r / 6)rdr 



+ 0(/3) 



-alog- + 0{P) 
log^ + 21og^ + 21n2 



+ 0(11 



(3.17) 



In the appendix some properties of ps are given when S is small, which will 
be useful for the next section. From expression (3.17), we obtain the following 
asymptotic expression for the mean and the variance of the number of bounds 



1— ^ ^^^^4+0(1) 



ms 



1 -ms 



mi 



(logi)(logi + 21ogi + 21og2) 



V 



(3.18) 



+ 0(1)(3.19) 



These expressions are valid for £ << 1 fixed, but are uniform in /? for /? << 1. 



4 Estimation of Mean First Passage Time 

E{T^\T^ < r|}. 

In this section, we give an asymptotic estimate of the mean time to hit 
denoted by T^, conditioned on the event that {T^ < Tj}. As described in 
[18], the conditional process of a Brownian motion conditioned on {T^ < Tj} 
satisfies the following stochastic differential equation 

dX*{t) = 2D^^4§^^dt + V2DdW. (4.1) 

On the outer boundary dD{R), the process X* is reflected (resp. absorbed) 
exactly where X is reflected (resp. absorbed). We denote t^{x) — E{Tx\Tx < 
T^}- It satisfies Dynkin's equation [10] which here is a degenerated elliptic 
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partial differential equation, with mixed boundary values: 

DpsAt"^ + 2DVt^ -Vps = -psonAs, (4.2) 

-^(x) — on dflr, 
on 

t^{x) = OondQa, 

where dQr, dfla are respectively the reflecting part and the absorbing part of 
the outer boundary. We remark that no boundary conditions are needed to be 
given in the inner circle dD{S) because Vp ■ n = ^ > and this is exactly the 
Fichera conditions, discussed in [22, 23], where boundary conditions cannot be 
given. 



4.1 Asymptotic expansion of the mean first passage time 

When the radius 6 of the inner circle is small, we obtain an explicit asymptotic 
expansion of the mean conditional time i^, solution of equation (4.2). We first 
derive an asymptotic solution by gluing two solutions: 1) when the initial point 
X is far from the inner-circle, the solution looks like the mean exit time when 
the drift term is set to zero. This solution is called the outer-solution and has 
been estimated in [24] with similar boundary conditions. 2) When the initial 
point X is chosen near the inner-circle, the solution can be approximated by 
a radial function. The approximation is valid in a boundary layer and has to 
match the radial part of the outer-solution at least . 



4.2 Outer-solution 

We now provide a construction of the outer-solution to equation (4.2). We 
start with the expansion of p, which depends on the parameter 5. 



ps{r, 9)^1- a{S, e)4>s{r, 9). 
By using the appendix, we obtain the following expression 



(4.3) 



n=l 



5(3 



cos(n^)+log(^)+41og(^))(4.4) 



Equation (4.2) can be written outside the boundary layer: Us — {r\rs < r < 
R}, 

D(l - a(5, e)(t)s{r, 9))At^ - 2a(5, e)V • M^, 0)^t^ = -(1 - a(5, e)(t>5{r, 9)) on 



dn 



t^{x) = on dQ,r: and t^{x) = on dQa- 
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We look for a regular asymptotic expansion of the solution 



t'^ix) = u{x) - a{5, e)ui{x) + 0{a\S, e)). (4.6) 

Using expression (4.6) and the behavior of p as 5 goes to zero, in the closed 
domain D{R) (see appendix), we obtain from equation (4.5) that u satisfies 

DAu = -1 on D(R), (4.7) 

— = on dflr, 
on 

u = on dQa- 

We can now use the result of [9, 24] to compute the leading order term of u. 
For X that does not belong to a boundary layer near the absorbing boundary 
dQa, the asymptotic expansion of u is given by 

u{x) = ^ ( In(^) + In 2) + 0{e) on Ug. (4.8) 

We conclude that u does not depend on the variable r and 6 at the first order 
in S and e. Because m is a solution of equation (4.7), outside a boundary layer 
of dQa: the derivatives |^ and ^ are small at the first order in S and e. We 
now consider the first order term in a in equation (4.5), which satisfies the 
equation 

DAui - 2V0Vu = 
^(x) = on dD.r, and ui{x) = on dfla- 

Outside the boundary layer near dQa, the leading order term is given by 
Ui{x) = o(l) (for the variable 6). We get that 

t"^{x) — u{x) + o{a{5,e)) on Us 

= ^(ln(-)+ln2)+o(a(5,e)) onC/5. (4.9) 



4.3 Asymptotic solution inside the boundary layer 

We now provide an asymptotic expansion of the mean time inside the bound- 
ary layer BLs — {S < r < r^}, where rs = S{1 — ^^^l^y)- The size of the bound- 
ary layer is given by formula (6.19) in the appendix. To estimate t^, we use 
the expression of the conditional probability ps- As described in the appendix 
(6.2.2), it can be approximated by a radial function, so that equation (4.2) 
becomes 

D({ty + -{t^y)+—{t^)'^ - -I on BLs, (4.10) 
V r J PS or 
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and the function r has to be glued continuously at r = r^, using that 

Ps{r, e) = _^iM (r-5)+ a{P, e)0 (r - 5) . (4.11) 
The leading order term is 

2D dps _ 2D 
Ps dr r — 5^ 

we have 



(4.12) 



D{{ty + -{t^y + 2^{t^y) = -ioubls. (4.13) 

r r — d 

The solution is given for r > (5 by 

('^)'(-) = 7(73^ (c-Ij(— -2*3 + 4') ("^) 

where the constant is chosen such that the function t"^ is integrable. Thus the 
numerator vanishes for r — 5. Actually r — 5 has to be a third order zero of 
the numerator and 

(t^)'(a.) = -^(r + <5/3). (4.15) 



Thus, 



where the constant C is determined by the matching condition: 



We get 



t^irs) = ^(In(^) +ln2) +0(a^(5,e)). (4.17) 



C^-m^{\n{^)+ln2^ (4.18) 



and for x e BLs, 

D V V"^e' ' '"~y ' 8q;2 8 ' 6 ' 12 



(4.19) 
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4.4 Asymptotic estimate of the average Mean time (r) 



We now compute asymptotically the mean time (r) defined in equation (2.10) 

by 

(r) = 7 y = AIM,, (4 20) 

where we recall that (r)^^ = ^^^^^^^ J^^t^{x)ps{x)dx. Vol{As) = 7r(i?^ — 6'^). 
To compute (t)^^, we decompose the domain As — BLg U {As — BLg). Using 
the previous computations for the outer solution (4.9) and the boundary layer 
solution (4.19), we compute each term separately and we get: 

kr)., = ln( j) + ln2) + + ln(-)tf + ^ !.(-) ln(- 

By taking into account the leading order term only, using the expressions of 
777.5 and p, we obtain after some computations that 

(r) ^ ^(ln(i) +ln2) - ^ ln(ln(i)) W(i)/3^ (4.21) 

We conclude that the boundary layer has very little influence on the leading 
order term at the order a. 



4.5 Computation of the mean time (T) for a molecule 
to hit the boundary dD{5) before exit 

We now provide an estimate of mean time it takes for a Brownian molecule 
to enter into the binding site domain D{5), before exit. To estimate the mean 
time (T), we observe that E{Tg\T§. < T^}qs{x) + E{T^\T^ < Ti}ps{x) is 
the mean time it takes for a Brownian molecule to exit the domain Q when it 
can either be absorbed in the inner disk or in the small absorbing boundary 
on the outer circle. If we denote the mean time by 

w{x) = E{Ti\Ti, < T^}qs{x) + E{T^\T^ < T^-jpsix) (4.22) 

then w satisfies the following PDE [10], 

in Q, 

D 

on dQa U dQs 
on dQr 



Aw = 

w — 
dw 

dn 



17 



We substitute w{r, 9) — v{r, 9) + ^^j^ , and solve the problem for v, solution 
of 



At; = 





in Q, 


V — 





on dO-a 


dv 
dn 


R 

2D 


on dVtr 


V — — 


AD 


on dVts 



We proceed as in the previous section by expanding v in Fourier Series 

CO „ 

^(^'^) = f + E K(^)" + ^"(7)"] '^'^''^^ + (^)- (^-2^^ 

n=l 

In a similar way as we estimated p and q (described in appendix A), from the 
boundary conditions on we get 

Y + E («"/^" + (;^) ") + 7log/3 = (4-24) 



from which we get 



ao = -^^^-27 log/? (4.25) 



With the definitions cq = Oq, c„ = a„(l — i/„ = ^^stt the boundary 

conditions on dVLridVLa have the form of double series equations 

00 

? + E TtV = °' ^ek-^,^] 

00 

^ nc„ cos(n6^) = (^^-7), 6'e[0,7r-£]. 

n=l 

The solution of these equations follows the steps described in the appendix 
and the asymptotic approximation of the coefficient cq is given by 

Co = 2 - 7) {2 log i + 2 log 2 + 0(/32)^ . (4.26) 

Using the last equation and (4.25) we obtain an equation for 7 

" 27 log /? = 2 f ^ - 7) (2 log - + 2 log 2 + 0(/?2, £)^ , 



2D ' 7 V e 
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which gives that 



and 



_ 41og^ + 41og2 + l + 0(/j^g) 
^ 2L>21og^+41ogi + 41og2 + 0(/52)' ^' ' 



^ (21ogl-l)(21ogl + 21n2) 

D 21ogi + 41ogi + 41og2 + 0(/32)- 



Remarks 



Averaging w with respect to a uniform distribution 

/ (21ogl - l)(logl + log2) - llogl - llog2 + llogl - i + 
"2^1^ logi + 21ogi + 21og2 + 0(/52) 

When /5 — > 0, we obtain the mean time to exit for only one small hole located 
on the boundary of the domain: this result agrees with the result obtained in 
([24], eq.(1.3)). In the limit S « e, we have the following approximation for 

M, 

(^)^:^((logi + log2) + i}, (4.29) 
while for £ << 5, we have 

M^^{(log^-^)-^}. (4.30) 
Final computations 

To compute E{td) it is sufficient to estimate the average q{x)E{T^\T§. < T^} 
defined by equation (2.11) 

Ja, <iiy)dy 

Starting from equation (4.22), we get 

q5{x)t^{x) = w{x) - p5{x)t^{x), 
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thus, 



J t^(x)q,{x)dx (u.) ms 



Using expression (3.17), we have 

21n(i) + 21n2 



1 /■ , , 21n(i) + 21n2 

^y^^,.(,M, = l-™.= >|A_^_^0(1). (4.33) 



Finally, the expression of (T) can be computed from (4.29) and (4.21) and in 
the limit S << e, the leading order term is 

To derive the final expression for the Dwell time, we gather the computational 
results: using equation (2.13), we have 

l-m.^21n(l) + 21n2 

ms In(^) 

Finally, we obtain from equation (2.12) that the dwell time in the approxima- 
tion S « s, 



(4.36) 

i?2 1 21n(i) + 21n2/ 1 i?2 l \ 

— (ln(-)+ln2)H + -{log-} +o(l). 

The expression of £'(rD) when e << 5 is more delicate and involves a different 
boundary layer analysis than the one given in the appendix. 



5 Discussion and conclusion 

In this article, we have computed the mean time spent by a molecule inside 
a microdomain Q, when it can interact with some binding sites, represented 
as a connected sub-domain D{5) C Q. This simplified assumption ignores the 
scattered distribution of the binding sites but leaves one free parameter 6. A 
rational way to choose the radius S is to equal the length 2^5 with the sum of 
the potential well boundary length of each binding site. 
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More generally, the problem of estimating the Dwell time of a molecule in 
fl when there are already many other independent molecules, is much more 
involved. The reason is that even if these molecules are not interacting, they 
are coupled through the competition for the binding sites. However, if the 
mean number of bound receptors is fixed and the variance is small enough, the 
Dwell time formula (4.36) can be apphed, by choosing for 5, a value that is 
related to the amount of free binding sites available (see below). In general, 
when the number of bounded molecules is fluctuating, a different model is 
needed to account for the fluctuations. In that different derivation of 

the Dwell time is needed. 

5.1 Receptor trafficking at synapses 

The explicit expression of the Dwell time can be used to describe receptor 
dynamics at synapses. A synapse is a micro-contact between two neurons, 
involved in signal transmission. In the past decades experimental observations 
have revealed that the molecular composition of the postsynaptic part of a 
synapse depends on the history of the neuronal activity [13, 12, 25, 14, 17, 15, 
26]. Although it is a difficult problem to predict the chemical organization of 
a synapse, it has been found experimentally that the type and the number of 
receptors do not appear randomly, but are well regulated during specific synap- 
tic plasticity protocols, such as Long Term Potentiation (LTP). During LTP 
[13], the number and the type of receptors can change, whereas extra- synaptic 
receptors diffuse inside a specific microdomain called the postsynaptic density 
(PSD). Recently, the concept has emerged that receptors are constantly mov- 
ing on the neuronal surface, in and out of the synaptic domain [15]. Moreover, 
receptors are also cycling between the cell surface and intracellular pools. Ac- 
cording to some recent experimental results [14], the movement of the receptors 
at synapses can be approximated by a random walk in the heterogeneous PSD. 
Moreover, it has been observed that receptors can become confined for random 
times [14] and can also be bound to scaffolding molecules. The PSD delimits 
a bounded domain fl containing small openings at the boundary, where recep- 
tors can be exchanged with the rest of the dendrite. The PSD contains many 
fundamental molecules, scaffolding proteins, receptors, kinases, and many oth- 
ers required for the normal synaptic activity. The PSD is thus a place where 
the synaptic molecular machinery is concentrated. Scaffolding molecules may 
be used to anchor receptors and/or to change the biophysical properties of 
receptors. Because scaffolding molecules are scattered inside the domain Q, in 
the present model we have replaced the complex organization of the PSD by a 
homogenized domain (Figure 1). We can apply the result of the first part of the 
present paper to obtain some estimates of the mean time spent by a receptor 
inside the PSD, when it can bind with scaffolding molecules. Estimate (4.36) 
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of the Dwell time E{td) depends on the size 5 of the scaffolding domain and is 
a good approximation when the amount of free scaffolding molecules is small 
5 «1 and when the corral barrier, measured by £, made by the impenetrable 
molecules around the PSD, is also small. 



5.2 Dwell time of a single receptor in a synapse con- 
taining many other receptors 

To estimate the mean time a receptor stays inside the PSD, we consider the 
case where the influx J of receptors entering the PSD is fixed. The Dwell time 
can be computed from formula (4.36) once the value of the radius 5 is known. 
5 is indeed the radius of the disk containing the free scaffolding molecules. To 
obtain an estimate of 5, we consider the dynamics where each injected receptor 
can escape with a rate r, approximated by formula (1.1). When the number 
of injected receptors is balanced by the exiting one from the PSD, the total 
number of free receptors inside the PSD is given by 

[R] = Jr. (5.1) 

(for a detailed analysis see [27]). Under the steady state assumption, the 
number of free scaffolding molecules can be estimated by using the law of 
chemical reaction (1.10) which says that 

_ [R-S] 

where [R — S] is the number of bounds. [R] (resp. [S]) the number of free 
receptors (resp. scaffolding molecules) K is the equilibrium reaction constant 
per unit area. If [Sq] denotes the initial number of scaffolding molecules, the 
conservation of matter (the total number of scaffolding molecules is conserved) 
and equation 5.2 implies that 

[S-o] = [5] + [i?-5] = [5] + i^[i?][5] (5.3) 
= {l + KJr)[S] (5.4) 

Thus, if [S] (resp. [Sq]) occupies a surface 7r5^, (resp. tt^q), then 

, . (5.5) 



5o can be related to the concentration cq of scaffolding molecules by cq — 



[So] 

To obtain the mean time estimate, expression (5.5) of 5 should be used in the 
Dwell time formula. 
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Expression (5.5) depends on the equilibrium constant K, which in fact 
depends imphcitly of d. Thus to provide a more accurate value of 6, we note 
that the forward binding rate ki represents the rate of arrival of a free receptor 
to a free binding site and can be approximated by 

h ~ (5.6) 

where (T) is given by equation (4.31). Using the previous considerations and 
the conservation law 5.3, we have 

[So] = |S| + (B-S) = lS] + A-lfl][S] (5.7) 

where r is approximated by expression (1.1) and (T) (4.34) depends on 5, 
which is denoted by {T){S). Using equation (5.5), 5 is solution of equation 

1 + 



In practice, the value of 5 obtained by solving equation (5.9) can now be used 
to estimate the Dwell time, the mean number of bounds and the other mean 
times. 

Finally, when the total number of receptors inside Vt is small, the radius 
5 cannot be approximated by a constant, rather it is fluctuating and in that 
case, the present approach needs to be adapted. A variant model of the PSD 
has been proposed in [27], based on a Markovian approach, which uses an in- 
and out-flux of receptors. At steady state, under some approximations, some 
estimates of the mean and the fluctuation of the total number of bound recep- 
tors are obtained, however the computations depend on an a priori expression 
of the forward binding rate. 



6 Appendix 

6.1 Solution of double series equation 

In this section we give the mathematical details needed to solve the double 
series equation (3.8)-(3.9). The method was already used in [24], [19]. Starting 
with the equations 

f + Er=i TT^ cos(ne) = 0, ^e[7r,7r-£] (6.1) 
« + Er=i^c„cos(ne) =0, ee[0,7r-£], (6.2) 
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we are interested in computing a and the coefficients {q}. Let us define the 
function h{9) for ^ e [0, tt — s] by 



^ ^ + ^ 7e -^cos(t^) — cos(t) 

Using this definition and the Fourier coefficient formula ([24]) we can write 

Cn h{t)[Pn{cOs{t)) + Pn-l{cOs{t))]dt, (6.4) 

Co = V2 / h{t)dt, (6.5) 

where Mchlcr's identity for Legendre polynomials gives ([19] ch.2) 

V2 r cos(n + |)d^ 
P„(cos(ii)) = — / =. 

^ Jo ■\/cos{9) - cos(w) 

Integrating equation (6.2), using the expression for the coefficients (6.4) into 
the integrated equation and the identity 

^ E~ :{/'»(cos(t)) + P„-:(coK*))) s^uO) = 

where H is the Heaviside function, we get the integral equation 

Jo ^/cos[t) - cos[9) Jo cos(|j 
where the kernel function K/^ is 

Kp{e,s) = Vg,(P,(cos(t) + P.(cos(t))). 

v2cos(|) ^ 

Using the first elements in the infinite series we can give an asymptotic expan- 
sion to 

Kp{t, s) = 2/3^ cos2(s) sin(i) + 0(/3^). (6.7) 
Equation (6.6) can be transformed according to the Abel transform 



Jo 

9{t) - 



g{t)dt 



{p{x) - pit)^ 
sin(7rg) d /"* p'{x)f{x)dx 



TT dtj^ {p{x)-p{t)y-<i' 
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where p{x) is monotonic increasing function. Transforming eq.(6.6), with q — \ 
and p{x) — — cos{x) we get 

h{Q) - / K^{e,t)h{t)dt = -- / , (6.8) 
Jo ^ Jo y cos t — cos 9) 

with ^^Cf S) - -i-^ K^K^)sin(»)) „ , n _ 1 d_ f-e 2atsm{^t)dt 

then equation (6.6) has the form (7 — X;3)/i = z. This is the Predholm integral 
equation, which can be approximated using the infinite sum 

h = z + Kpz + Kjz + . . . (6.9) 

Computing the first elements in the sum gives an approximation to h{9) which 
in turn, by substituting the approximation into the expressions for c„ gives an 
approximation for the coefficients. To calculate cq up to 0{e,l3^) we have to 
evaluate the integral 

j z{t)dt + j j K{s,t)f^z{t)dtdsj. 

The computation is made in [24] and gives 

Co 2a;(2 log - + 2 log 2 + 4/3^) . 
e 



6.2 Properties of the conditional probability ps 



6.2.1 Inside As 



qs is solution of equation (3.1) and depends on the parameter S. 



Proposition 1 The sequence ps converges uniformly to the constant 1, on any 
compact set K strictly contained in the domain A^. More specifically, for any 
K <Z As, there exists a constant C > such that for all r,9 & K, 

\psir,9)-l\ < CaiP,e) 



\^Psir,9)\ < Ca{P,e) 
where a{(3,e) is defined by (3.12) and 



<«(/?,£) = a(-^,£)<^ (6.10) 



which tends to zero when 5 goes to zero. 
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Proof. To estimate ps we consider expression (3.2), which can be written as 



n=l 



2n 



f \ n 



cos{n9) + a{P, e) log j . 



where a{P,e) is given by expression (3.12) and we recall that 



1 + H„ 



Ct7, 



V2 



/ h{t)Pn{cOs{t)) + Pn-l{cOs{t))]dt. (6.11) 

Jo 



By definition Hn = 0{j3^^) and P„ is the Legendre polynomial. For x G [—1, 1] 
and all n, |Pn(2^)| ^ 1- Moreover, using formula (6.9), we obtain that 



h{t) = z{t) + 0{P^) = - 



2a d 



usmu 



TT dt Jq -^/coslf^^cosT 



du. 



(6.12) 



Thus we conclude from the exphcit formula (6.11) that there exists a constant 
C > such that for all n 



\cn\ < Ca{(3,e) 



(6.13) 



which is independent of n. We denote c„ = c„q;(/3,£). We can now obtain the 
desired estimates. For tq < r < Ri < R, we get. 



l-p5(r,^) = l + a(/3,£)(5^- 
Thus, 



2n 



where 



\l-ps{r,9)\<I + II, 



I<\l-a((3,s)\og (^)l<^ 



,l n(i?i/ro) + 0(l) 
\n{R/5) + 0(1) 



(6.15) 

(6.16) 
(6.17) 



The last part of these inequalities uses the asymptotic expression of a given 
by formula (3.12). These inequalities show that on any compact set of As, 
1 — Ps converges uniformly to zero at a rate j^y/^J- Similarly, to estimate the 
gradient we estimate separately ||^| and |^|^|- First, 



dps 
dr 



.«)-(^.<f:T^%.[(9^)-(^)]»^(««)^i). 
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we get 



dr 

and second, 
r de ^ 

and 



n=l 



n=l 



r 9^ 



oo 



n=l 



5^ 



6.2.2 The boundary layer 

We consider now the behavior of ps near the the inner circle r = d, where ps 
vanishes. We expect to see a boundary layer near the circle r = d, as d goes to 
zero. Ps is a regular function and a Taylor expansion in the r- variable inside 
expression (6.14) gives that 

a{P, e 



n=l 



r — 5 



The leading order term is radial and we get 



P5{r, 9) = (r - 5) + e)0 (r-d). 



(6.18) 



To estimate the size of the boundary layer, we look at the value rs such that 
Ps is of the order one. Using the previous approximation of ps (6.18) we get 



rs ^ 5 



5 + 5(lni + 21n- + 21n2 + lni? + 0(/32)). 
6 e 



(6.19) 
(6.20) 



6.2.3 Convergence of ps in the neighborhood of dD{R) 

The sequence of harmonic function ps converges uniformly on any compact set 
of D{R) — {0} to the function 1. We show now that this is also the case in 
the neighborhood of dD{R). We use the elliptic estimates given in [28] derived 
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at the boundary. We can also directly extend as a harmonic function ps to a 
neighborhood T of D{R). We have that 



sup iVp^l < C (6.21) 
{xeK} 

for any compact K in T. By Ascoh's theorem, there exists a subsequence of 
Ps converging in C°'". We conclude since ps converges to 1 on any compact set 
of D{R) — {0} that the sequence ps converges to 1 on T — {0} and is a weak 
solution of 

on D(R) - (6.22) 

1 on dQa: 

on dD{5), 

Because < m < 1, the solution can be extended up to and the only solu- 
tion of equation (6.22) is the constant 1. We conclude that the sequence ps 
converges to the constant one uniformly on D[R). 



Au = 
u{x) — 

> - 
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